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Abstract 

Within a theoretical study using both HYLANDER and RQMD, we revisit Bose-Einstein 
correlation measurements of 200AGeV S+S and 160AGeV Pb+Pb. Transverse flow does not 
show up in the Bose-Einstein measurements. The decrease of the effective transverse radius 
parameters with increasing transverse average momentum of the particle pair is largely due to 
the effect of resonance decays. 



*E. Mail: schlei@t2.LANL.gov 
t E. Mail: nu_xu@lanl.gov 



1 



In the ongoing search for the quark-gluon plasma (QGP), Bose-Einstein correlations (BEC) 
of identical bosons have become an issue of great current interest. BEC serve as a tool in the 
determination of radii (spatial extensions) and lifetimes (temporal extensions) of hadron emission 
zones generated in relativistic nucleus-nucleus collisions However, the inverse widths of BEC 
functions do not have a simple interpretation in terms of those radii and lifetimes. The inverse 
widths or "correlation radii" have to be considered as combinations of space-time moments from 
space-time densities (c/., e.g., ref. Q and references therein). The phase-space distribution, or the 
source function g(x,p), does not only depend on the space-time point x and the four-momentum 
p of the emitted particles, but also on the temporal evolution of the hadron source (fireball) and, 
more importantly, on the correlations among x and p Furthermore, it is known that BEC are 
sensitive to the effects of the decay of unstable particles (resonances). Experimentally obtained 
BEC functions are always averaged over certain phase-space regions and are therefore in general 
sensitive to the specific kinematical regions under consideration. Because of the complications 
involved in BEC, it is necessary to account for all known effects before we can draw conclusions 
about new physics which might emerge in relativistic nucleus- nucleus collisions. 

Recent theoretical studies [||-[[^] claim that the inverse widths of BEC functions in relativistic 
heavy-ion collisions at CERN/SPS energies near 200 GeV per nucleon show a mj_-scaling behavior, 
namely that the size parameters are inversely proportional to the square-root of the transverse 
mass, yjm±. This idea was proposed earlier by Sinyukov et al. || in 1987, in which the authors 
suggested that the longitudinal size parameter, R» , should follow the mj_-scaling: 



where Tj and yx are the particle freeze-out temperature and rapidity of the particle pair, 



respectively. The pair transverse mass is defined as m± = J m 2 + Kj_ with the average transverse 
momentum K± = (pi± + P2±)/2 of the particle pair. This relationship has been indeed observed 
in the longitudinal size parameter Rn from the sulphur-induced high energy collisions |J and the 
effect is attributed to the collective expansion. Neglecting the resonance decays, the authors in refs. 

suggested that even in the transverse directions, size parameters R ou t and R s ide also follow 
the same m^-scaling law. On the other hand, it has been shown that the effect of the resonance 



decay strongly affect the inverse widths of the BEC [Kj]-|l2]. 



It is the purpose of this short note to study the m_i_-scaling behaviour by using two realis- 
tic source models: HYLANDER [13] and RQMD |l4j| . This study includes effects of transverse 
expansion and decay of resonances on BEC of identical particles. We shall first compare HYLAN- 
DER results with recent measurements of S+S and Pb+Pb collisions. Then we will focus on the 
mj_-scaling issues by inspecting the mj_-dependence of the size parameters calculated from both 
models. 

The HYLANDER model belongs to the class of models applying (3+l)-dimensional relativistic 
one-fluid-dynamics (for other hydro dynamical models cf. [^] and refs. therein). HYLANDER pro- 
vides fully three-dimensional numerical solutions of the hydrodynamical relativistic Euler-equations 
Jig]. The model has been successfully applied to several different heavy-ion collisions at SPS en- 



ergies [17]. Once the initial conditions and the equation of state are specified, one obtains an 
unambiguous solution from the hydrodynamical equations. The choice of a freeze-out condition 
such as, e.g., that the expanding fluid reaches a fixed freeze-out temperature, Tf, determines the 
final space-time geometry of the hydrodynamically expanding fireball. In particular, for 200AGeV 
S+S and 160AGeV Pb+Pb central collisions, the calculation of single inclusive spectra of negative 
hadrons and protons and BEC of negative pions and kaons was performed using the formalism 



outlined in refs. |L(], [IT], 18]. Resonance decays were included in these calculations. The results for 



BEC of negative pions and kaons have been published in refs. [10|-[12],|15]. 

In order to extract effective hadron source radii, the results of the BEC calculations have been 
fitted to the Gaussian formQ which has been widely used by experimentalists for the presentation 
of their BEC data (for the choice of the variables, cf. ref. |20|j): 



Ciipulk) = 1 + A exp 



~ (q\Ri(K) + q 2 ide R 2 ide (K) + q 2 out R 2 out (K) 



(2) 



2 

In eq. (|2[) pi, p2 are the momenta of the two emitted identical particles, whereas K denotes the 

average momentum of the particle pair. The qi (i = ||, side, out) are components of the momentum 

difference vector q of the particle pair. It should be emphasized that in the present model A = \{K) 

does not represent the effect of coherence, but the momentum-dependent effective reduction of 

the intercept due to the contributions from the decay of long-lived resonances (such as ry, rf ... 

1 If one performs a fit to a BEC function in more than one dimension, eq. (^) does not represent the most general 
expression, because of the existence of an "out- longitudinal" cross term [^9) . 
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etc.) [[101 - 1 12 ],| 151- A characteristic property of particle production from an expanding source is 
a correlation between the space-time point where a particle is emitted and its energy-momentum 
[||. As a consequence, the inverse widths Ri(K) (i = \\, side, out) extracted from Bose-Einstein 
correlation functions show a characteristic dependence of the average momentum of the pair, K, 
and are therefore dependent on detector acceptances. 

Fig. 1 (solid lines) shows the calculations for the effective radii Rn , R s id e and R ou t as functions of 
the transverse momentum K± of the pion pair compared to the corresponding NA35 and preliminary 
NA49 data || respectively. In order to make a comparison of the calculated effective radii with 
the experimentally obtained ones, detector acceptances have been accounted for. In the case of 
200AGeV S+S the effective radii have been calculated as functions of K± at the rapidity of the 
pair dk = 4.0 — y cm ~ 1.0 . In the case of 160AGeV Pb+Pb the effective radii were calculated as 
functions of K± at yx = 4.5 — y cm ~ 1.6. The effective longitudinal radii Rn are evaluated in the 
longitudinal comoving system (LCMS), with 711 = cosh(y^). All of these calculations, which in the 
case of S+S had been true predictions, agree surprisingly well with the data. 

The effective radii of only the directly emitted (thermal) negative pions, ir~ , (dashed lines) are 
also plotted in Fig. 1. We can see that the effective radii of the directly emitted ir~ are smaller 
compared to those of all ir~ , i.e., including the pions originating from resonance decays, and follow 
a different i^-dependence. The effective radii for all pions are larger compared to the effective 
radii of directly emitted pions, because the unstable particles (resonances) have a finite lifetime. 
Within their lifetime the resonances can propagate to a distant location from their points of origin 
(the fireball). While decaying the resonances generate a cloud of pion emission locations around 
the actual fireball and therefore increase the effective size of the system. 

Fig. 1 shows that the -fTj_-dependence of the effective transverse radii of thermal pions is very 
weak while the effective transverse radii of all pions significantly decrease as a function of K±. 
In refs. [l5| it was explained that the hydrodynamical solutions express transverse expansion. 
The two systems each have a maximum value for the transverse velocities: for 200AGeV S+S 
we have a maximum transverse velocity u™ ax (S) = 0.43 whereas for 160AGeV Pb+Pb we have 
u™ ax (Pb) = 0.61. The maximal values for the transverse velocities correspond to fluid cells which 
have space-time rapidities 77 ~ 0.5 • ln(t + z/t — z) equal to zero (t and z are the time and the 
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longitudinal spacial coordinate, respectively). On the contrary, fluid cells with large space-time 
rapidity have very small corresponding values for the transverse velocities at freeze-out (cf. Fig. 6 
in ref. |l5[ ). We have checked that in the phase-space regions, where the BEC data of the NA35 
and NA49 Collaborations have been taken, the contributions to the BEC functions come mainly 
from fluid cells with almost vanishing transverse velocity at freeze-out. 

Transverse flow has two major effects on effective BEC radii. First there is an apparent re- 
duction of the effective longitudinal radii, because the fluid also expands transversely. Secondly, 
the transverse effective radii show a -fTj_-dependence. For the particular hydrodynamical solutions 
under consideration transverse flow results in a decrease of the transverse effective radii due to the 
effects of relatively fast transverse inwardly moving rarefaction waves. In agreement with ref. [18], 
the effective transverse radii of directly emitted negative pions plotted in Fig.l show almost no 
-Kj_-dependence. But the effective radii of all tt~ have a -/^-dependence which would indicate a 
present transverse flow in the data. Since there is no transverse flow in the considered phase-space 
region as we have argued above, the decrease of the effective radii for all 7r~ must have a differ- 
ent reason. In fact, for smaller transverse momenta we have relatively larger contributions from 
resonance decay to the BEC than at higher transverse momenta K±. Since the relative weight 
of pions being produced directly compared to those originating from resonance decays enters in 
the composition of BEC functions, the decrease of the effective radii of all pions with increasing 
K± can be understood in terms of the chemical composition at freeze-out as a function of particle 
momenta. Resonance decays produce a pionic cloud that surrounds the fireball. This pion cloud 
has in general larger spatial extensions (and a longer lifetime) than the fireball which emits the 
direct pions. Thus relatively larger resonance contributions result in relatively larger correlation 
radii and therefore, the decay of resonances is responsible for the decrease of the transverse radii 
with increasing K±. 

Let us now discuss mj_-scaling according to refs. If mj_-scaling would be present, we 

would get 

l\R\{K) ~ R 2 slde (K) ~ R 2 out (K) oc ^L. (3) 
The HYLANDER BEC calculations for the transverse radii R ou t and R s ide show an even stronger 



decrease with K± at midrapidity, i.e., yx = 0.0 (cf. Fig. 2 in ref. [1C] and Fig. 3a in ref. [12]). 
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In order to get a maximum effect due to a possible mj_-scaling we show in Fig. 2 effective radii 



calculated at yx = 0.0 and multiplied with Jm_\_/Tf as a function of the pair transverse mass m±. 
A present mj_-scaling would result in a constant behaviour of the plotted quantities, i.e., there 
would be no change with changing m±. Obviously, the effective radii multiplied with the factor 
m±/Tf do not show an independence of m±. Consequently, the effective radii R\\, R s ide and Rout 



(and those of only directly emitted tt~) are not consistent with a simple mj_-scaling behaviour in 
this hydrodynamical model. 

Up to now, we have based all of our arguments on the hydrodynamical model HYLANDER and 
simple mj_-scaling is found neither in the longitudinal direction, R», nor in the transverse directions, 
Rout and R S ide- In addition, we found that resonance decays increase the size parameters at low 
transverse momenta dramatically. 

To further understand the physics of the transverse momentum dependence of the size param- 
eters, we made an additional test with the widely used transport model RQMD (v2.2) [22[ |. Using 
the freeze-out phase-space distribution, g(x,p), we calculated the averaged pion size parameters 
R\\, Rout, and Rside as a function of the pair transverse mass m± for 160AGeV Pb+Pb central col- 



lisions, multiplied with ^Jm±/Tf and again at midrapidity, yx = 0.0. The results are shown in Fig. 
3 where the solid lines represent all pions and dotted lines are for rescattered pions. Rescattered 
pions do not originate from resonance decay at the moment of their last interaction. A temperature 
parameter of Tf = 140MeV is used here p2[ . The rescattered pions are similar to the thermal 
pions from the hydrodynamical model calculations (see dashed lines in Figs. 1,2), because they are 
emitted earlier. 

It is well known by now that, due to rescatterings in the heavy- ion collisions, collective flow 



manifests itself in the cascade type of calculations [23, 24]. In the Pb+Pb central collisions, the 
averaged transverse flow velocity at midrapidity is about 40% of the speed of light. If the m±- 



scaling [Q] is correct, one would expect a flat distribution of the scaled size parameters y m±/TfRi 
(i =[[, out, side) with respect the pair transverse mass, m±. However, Fig. 3 does not show such 
a strong flow effect as proposed in ref. ||]. It is worth noting, that in a recent publication p5[ on 
S+Pb central collisions one indeed observed a mj_-scaling. The corresponding RQMD predictions 
which accounted for the exact experimental acceptance were consistent with the data. However, 
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the experimental acceptance has a rapidity span larger than one unit and, at each transverse 
momentum, the corresponding rapidity is different. Such an experimental effect may explain the 
different m^-behaviour between the experimental results of NA35 || and NA44 pjj] and also the 
difference between our results in Fig. 3 and NA44. 

To summarize, we have revisited Bose-Einstein correlation measurements for 200AGeV S+S and 
160AGeV Pb+Pb. For the theoretical analysis we have used a hydrodynamical model (HYLAN- 
DER) and a cascade model (RQMD). Within the experimental acceptance, the hydrodynamical 
model results are consistent with the measured size parameters for both 200 AGeV S+S and 160 
AGeV Pb+Pb central collisions. However, at midrapidity, no mj_-scaling is seen in the results from 
both model calculations. Resonance decays seem to play an important role in the Bose-Einstein 



correlation radii. Our conclusion is also consistent with recent publications [26, 27 1. 
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Figure Captions 



Fig. 1 Effective radii extracted from Bose-Einstein correlation functions of identical pions as a 
function of the transverse average momentum K± of the pion pair at rapidity yx = 4.0— y cm ~ 
1.0 for 200AGeV S+S and at rapidity y K = 4.5-y cm « 1.6 for 160AGeV Pb+Pb, respectively. 
Solid lines correspond to all ir~ , i.e., including those from resonance decay, whereas the dashed 
lines correspond to directly emitted (thermal) ir~ . The theoretical results are compared to 
NA35 data [|, g| and preliminary NA49 data [§. 

Fig. 2 Effective radii extracted from Bose-Einstein correlation functions of identical pions mul- 
tiplied with ^Jrn±jTf as a function of the transverse mass m± of the pion pair at rapidity 
yx = 0.0 for 200AGeV S+S and for 160AGeV Pb+Pb, respectively. Solid lines correspond 
to all ir~ , i.e., including those from resonance decay, whereas the dashed lines correspond to 
directly emitted (thermal) ir~ only. 

Fig. 3 Scaled radii from the transport model RQMD (v2.2). A freeze-out temperature of 140 MeV 
is used in the figure. Similar to Fig. 2, one does not find the mj_-scaling law in these radius 
parameters. 
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Figure 1 
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